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ABSTRACT 

The most important evidence of non-thermal phenomena in galaxy clusters comes 
from Giant Radio Halos (GRHs), spectacular synchrotron radio sources extended over 
^Mpc scales, detected in the central regions of a growing number of massive galaxy 
clusters. A promising possibility to explain these sources is given by in situ stochastic 
reacceleration of relativistic electrons by turbulence generated in the cluster volume 
during merger events. Cassano & Brunetti (2005) have recently shown that the ex- 
pected fraction of clusters with radio halos and the increase of such a fraction with 
cluster mass can be reconciled with present observations provided that a fraction of 
20-30 % of the turbulence in clusters is in the form of compressible modes. 

In this work we extend the above mentioned analysis, by including a scaling of 
the magnetic field strength with cluster mass. We show that, in the framework of the 
reacceleration model, the observed correlations between the synchrotron radio power 
of a sample of 17 GRHs and the X-ray properties of the hosting clusters are consistent 
with, and actually predicted by a magnetic field dependence on the virial mass of 
the form B (x M^, with 6>0.5 and typical /iG strengths of the average B intensity. 
The occurrence of GRHs as a function of both cluster mass and redshift is obtained: 
the evolution of such a probability depends on the interplay between synchrotron and 
inverse Compton losses in the emitting volume, and it is maximized in clusters for 
which the two losses are comparable. 

The most relevant findings are that the predicted luminosity functions of GRHs are 
peaked around a power Piaghz ^ lO^'' W/Hz , and severely cut-off at low radio 
powers due to the decrease of the electron reacceleration in smaller galaxy clusters, 
and that the occurrence of GRHs at 1.4 GHz beyond a redshift z ^ 0.7 appears to 
be negligible. As a related check we also show that the predicted integral radio source 
counts within a limited volume (z ^ 0.2) are consistent with present observational 
constraints. Extending the source counts beyond z=0.2 we estimate that the total 
number of GRHs to be discovered at ~ mJy radio fluxes could be ^ 100 at 1.4 GHz. 
Finally, the occurrence of GRHs and their number counts at 150 MHz are estimated 
in view of the fortcoming operation of low frequency observatories (LOFAR, LWA) 
and compared with those at higher radio frequencies. 

Key words: particle acceleration - turbulence - radiation mechanisms: non-thermal 
- galaxies: clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

The intracluster medium (ICM) is a mixture of thermal and 
non-thermal components and a precise physical description 
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of the ICM also requires adequate knowledge of the role of 
non-thermal components. 

The most detailed evidence for non-thermal phenomena 
comes from the radio observations. A number of clusters 
of galaxies are known to contain wide diffuse synchrotron 
sources (radio halos and relics) which have no obvious con- 
nection with the individual cluster galaxies, but are rather 
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associated to the ICM (e.g., Giovannini & Feretti 2000; 
Kempner & Sarazin 2001; see Giovannini & Feretti 2002 for 
a review) . The synchrotron emission of such sources requires 
a population of GeV relativistic electrons (and/or positrons) 
and cluster magnetic fields on fiG levels. Evidence for rela- 
tivistic electrons (and positrons) in the ICM may also come 
from the detection of hard X-ray (HXR) excess emission in 
the case of a few galaxy clusters (e.g., Rephaeli & Gruber 
2003, Fusco-Femiano et al. 2004), and possibly from extreme 
ultra-violet (EUV) excess emission (e.g., Kaastra et al. 2003; 
Bowyer et al. 2004) . It is also believed that the amount of the 
energy budget of high energy protons in the ICM might be 
significant, due to the confinement of cosmic rays over cos- 
mological time scales (Volk et al. 1996; Berezinsky, Blasi & 
Ptuskin 1997; Enfilin et al. 1997). Nevertheless, the gamma 
radiation that would allow to infer the fraction of relativistic 
hadrons in clusters has not been detected as yet (Reimer et 
al., 2003, see Pfrommer & Enfilin 2004 for upper limit on 
this fraction). 

Shock waves are unavoidably formed during merger 
events; they may efficiently accelerate relativistic particles 
contributing to the injection of relativistic hadrons and 
of relativistic emitting electrons in the ICM (e.g., Ryu et 
al. 2003, Gabici & Blasi 2003). However the accelerated elec- 
trons have a short pathlength due to IC losses and thus they 
can travel a short distance away from the shock front, emit- 
ting synchrotron radiation concentrated around the shock 
rim (e.g., Miniati et al. 2001). Radio Relics, which are po- 
larized and elongated radio sources located in the cluster 
peripheral regions, may indeed be associated to these shock 
waves, as a result of Fermi-I diffusive shock acceleration of 
ICM electrons (Enfilin et al. 1998; Roettiger et al. 1999), 
or of adiabatic energization of relativistic electrons confined 
in fossil radio plasma, released in the past by active radio 
galaxies (Enfilin & Gopal- Krishna 2001; Hoeft et al. 2004). 

The most spectacular evidence of diffuse synchrotron 
emission in galaxy clusters is that associated to giant radio 
halos, Mpc-size radio sources which permeate the cluster 
volume similarly to the X-ray emitting gas. In this respect, 
two main possibilities have been investigated in some de- 
tail to explain that GeV electrons (and/or positrons) are 
present and able to radiate on distance scales larger than 
their typical loss lengths: i) the so-called reacceleration mod- 
els, whereby relativistic electrons (and positrons) injected in 
the ICM by a variety of processes active during the life of 
galaxy clusters are continuously re-energized in situ during 
the life-time of the observed radio halos (which is estimated 
to be ~ 1 Gyr, Kuo et al. 2004) and ii) the secondary elec- 
tron models, whereby electrons are secondary products of 
the hadronic interactions of cosmic rays with the intraclus- 
ter medium, as first proposed by Dennison (1980). Although 
the origin of the emitting particles in radio halos is still a 
matter of debate (e.g., Enfilin 2004), the above two models 
for the production of the radiating electrons (and positrons) 
have a substantial predictive power, which can be used to 
discriminate among such models by comparing their pre- 
dictions with observations. Although future observations re- 
main crucial to achieve a firm conclusion, at least as far 
as the few well studied clusters and the analysis of statisti- 
cal samples are concerned, present data seem to suggest the 
presence of in situ particle-reacceleration mechanisms ac- 



tive in the ICM (e.g., Brunetti 2003,04; Blasi 2004; Feretti 
et al. 2004; Hwang 2004; Reimer et al. 2004). 

Radio observations of galaxy clusters indicate that the 
detection rate of radio halos shows an abrupt increase with 
increasing the X-ray luminosity of the host clusters. In par- 
ticular, about 30-35% of the galaxy clusters with X-ray lu- 
minosity larger than 10'*^ erg/s show diffuse non-thermal 
radio emission (Giovannini & Feretti 2002); these clusters 
have also high temperature (kT > 7 keV) and large mass 
(> 2x 10^^ Mq). Furthermore, giant radio halos are always 
found in merging clusters (e.g., Buote 2001; Schuecker et al 
2001). Although the physics of particle acceleration due to 
turbulence generated in merging clusters has been investi- 
gated in some detail (e.g., Schlickeiser et al. 1987; Petrosian 
2001; Fujita et al 2003; Brunetti et al. 2004; Brunetti & 
Blasi 2005) and the model expectations seem to reproduce 
the observed radio features and possibly also the hard X- 
rays (e.g., Brunetti et al. 2001; Kuo et al. 2003; Brunetti 
2004; Hwang 2004), a theoretical investigation of the statis- 
tical properties of the Mpc diffuse emission in galaxy clusters 
in the framework of these models has not been carried out 
extensively as yet. In particular, the fact that giant radio 
halos are always associated to massive galaxy clusters and 
the presence of a trend between their radio power and the 
mass (temperature. X-ray luminosity) of the parent clusters 
may be powerful tools to test and constrain present models. 

In a recent paper Cassano & Brunetti (2005; hereafter 
CB05) have modelled the statistical properties of giant radio 
halos in the framework of the merger-induced in situ particle 
acceleration scenario. By adopting the semi-analytic Press 
& Schechter (1974; PS74) theory to follow the cosmic evolu- 
tion and formation of a large synthetic population of galaxy 
clusters, it was assumed that the energy injected in the form 
of magnetosonic waves during merging events in clusters is 
a fraction, rjt, of the PdV work done by the infalling sub- 
clusters in passing through the most massive one. Then the 
processes of stochastic acceleration of the relativistic elec- 
trons by these waves, and the ensuing synchrotron emission 
properties, have been worked out under the assumption that 
the magnetic field intensities, ICM temperatures and parti- 
cle densities (both thermal and non-thermal) have constant 
volume averaged values (within 1 Mpc^). The main find- 
ings of these calculations is that giant radio halos are nat- 
urally expected only in the more massive clusters, and the 
expected fraction of clusters with radio halos (at redshifts 
z< 0.2) can be reconciled with the observed one under viable 
assumptions (?7t — 0.24 — 0.34 ). Specifically, the probability 
to form giant radio halos in the synthetic cluster population 
was found to be of order 20-30 % in the more massive galaxy 
clusters (M > 2 x 10^^ Mq), 2-5 % in M ~ lO^'^ Mq clusters, 
and negligible in less massive systems. Such increase of the 
probability with the cluster mass is essentially due to the 
increase of both the energy density of turbulence and of the 
turbulence injection volume with cluster mass (see CB05). 

The present paper is a natural extension of the CB05 
work, the most important difference being that here we 
adopt a scaling law between the rms magnetic field strength 
(averaged in the synchrotron emitting volume) and the virial 
mass of the parent clusters, B oc M'' . We carry out a de- 
tailed comparison between statistical data of giant radio ha- 
los currently available and model expectations as derived by 
adopting the CB05 procedures. 
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Figure 1. Correlation between the radio power at 1.4 GHz and 
the X-ray luminosity between [0.1-2.4] kev for the GRHs. 



In Sec. 2 we collect radio and X-ray data for well known 
giant radio halos from the literature and derive radio-X-ray 
correlations. In Sec. 3 we investigate the possibility to match 
the observed radio-X-ray correlations for giant radio halos 
with electron acceleration models. This comparison provides 
stringent constraints on the physical parameters in the ICM, 
in particular for the magnetic field in galaxy clusters. In 
Sec. 4 we derive the expected probability to form giant radio 
halos as a function of M„ and z. This is done by adopting 
the same values of the physical parameters which allows to 
account for the observed radio-X-ray correlations. In Sees. 5- 
6 we finally calculate the expected luminosity functions and 
number counts of giant radio halos. 

As in CB05, we focus our attention on giant radio ha- 
los only (linear size ~1 h'^g Mpc, GRHs elsewhere). The 
adopted cosmology is: ACDM {Ho = 70 Km s"^ Mpc~'^, 
no,m = 0.3, Qa = 0.7, 0-8 = 0.9). 



2 OBSERVED CORRELATIONS 

In this section we discuss the observed correlations between 
the X-ray and the radio properties of clusters hosting GRHs. 
We collect galaxy clusters with known GRHs from the lit- 
erature obtaining a total sample of 17 clusters. In Tab. 
we report the radio and X-ray properties of this sample in 
a ACDM cosmology. In order to have the best estimate of 
the X-ray temperatures we select results from XMM-Newton 
observations when available, otherwise we use ASCA results 
or combine ASCA and Chandra information. We investigate 
the correlations between the X-ray and the radio properties 
of the selected clusters by making use of a linear regression 
fit in log-log space following the procedures given in Akritas 
& Bershady (1996). This method allows for intrinsic scatter 
and errors in both variables. 



Figure 2. Halo radio brightness at 1.4 GHz, normalized to the 
radio brightness of A3562 (which is just visible in the NVSS), 
versus X-ray luminosity between [0.1-2.4] keV. Different symbols 
indicate GRHs (filled circles) and smaller radio halos (open cir- 
cles) visible in the NVSS. Asterisks mark A2256, which falls below 
the NVSS surface brightness limit, and 1E50657-558, which falls 
below the declination range of the NVSS. Large circles mark ob- 
jects visible in the NVSS and in the redshift range z ~ 0.15 — 0.3. 



2.1 Radio Power— X-ray luminosity correlation 

The presence of a correlation between the radio powers and 
the X-ray luminosities is well known (Liang et al. 2000; Fer- 
etti 2000, 2003; Enfilin and Rottgering 2002). 

In Fig0 we report the correlation between the X-ray 
luminosity (in the 0.1-2.4 keV energy band) and the radio 
power at 1.4 GHz (A. 4) for our sample of GRHs. The fit 
has been performed by using the form: 



log 



Pi. 



3.16- 1024/1; 



= Af+bf lo] 



L. 



1045 /l"' ^ 



(1) 



0.159 ±0.060 and 



where the best fit parameters are: Af 
bf = 1.97 ±0.25. 

Our findings are consistent with those of Enfilin and 
Rottgering (2002) who used 14 clusters with radio halos 
and found a correlation of the form Piaghz oc I/jf^**. Us- 
ing 16 clusters with GRHs Feretti (2003) found a correla- 
tion between the X-ray bolometric luminosity and the ra- 
dio power at 1.4 GHz of the form Piagh. oc (l6^')1-8±0-2. 
A consistent result is obtained with the data in Tab. 1 

[PlAGHz OC [Lx ) )■ 

Although the trend in FigQ appears quite stringent, 
one may wonder wether it could be affected by observational 
biases. It should be pointed out that most GRHs have been 
discovered by follow ups of radio halo candidates mostly 
identified from the NRAO VLA Sky Survey (NVSS, Condon 
et al. 1998). 

In Fig|5] we plot the average radio surface brightness of 
GRHs at 1.4 GHz, normalized to the radio brightness of 
A3562 (open circle with the smallest X-ray luminosity), 
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Table 1. Radio and X-ray properties of cluster with GRHs ( linear size ^ 1 hrg Mpc) in a ACDM cosmology. In Col.(l): Cluster name. 
Col. (2): Cluster redsliift. Col. (3): Cluster temperature given in keV. Col. (4): X-ray luminosity in the energy range [0.1 — 2.4] keV in unit 
of /lyQ 10'*'* erg/s. Col. (5); Bolometric X-ray luminosity in the energy range [0.01 — 40] keV in unit of h^Q 10*** erg/s. Col. (6): Radio 
power at 1.4 GHz in unit of /iZ, 10'^'* Watt/Hz. Col. (7): Large Linear Size (LLS) of the Radio Halo is in /ir„ kpc. Ref. for the temperature 
data in brackets: (Z04) Zhang al. 2004 (XMM); (WOO) White et al. 2000 (ASCA); (M96) Markevitch 1996 (ASCA); (m) mean value 
between Mushotzky & Scharf 1997 (ASCA) and Govoni et al. 2004 (Chandra); (e) Ebeling et al. 1996 (from L^-T relation) ; (D93) David 
et al. 1993 (Einstein MPC+ Exosat -|- Ginga); (M98) Markevitch et al. 1998 (ASCA); (ml) mean value between Z04 and Pierre et al. 
1999 (ASCA data); (H93) Hughes et al. 1993 (GINGA). Ref. for the X-ray luminosities in brackets: (B04) Boehringer et al 2004, (E98) 
Ebeling et al 1998, (E96) Ebeling et al 1996, (T96) Tsuru et al 1996, Ref. for the radio data in brackets: (LOO) Liang et al. 2000 (ATCA) 
(FOO) Fcrctti 2000, (B03) Bacchi ct al 2003, (GFOO) Giovannini & Feretti 2000, (V03) Venturi et al 2003, (GFGOl) Govoni ct al. 2001, 
(G05) Govoni et al. 2005, (FF03) Fcrctti et al. 2001, (m2) mean value between Kim ct al. 1990 and Dciss ct al. 1997 



cluster's 


z 


T 


Lx 


Lboi 


Pi A 


LLS 


name 




[keV] 


[10*4 erg/s ] 


[1044 crg/s ] 


[102'4 Watt/Hz] 


[Mpc h-,l] 


1E50657-558 


0.2994 


13-59lo58(Z04) 


23.322 ± 1.84(B04) 


88.619 ±7.00 


28.21 ± 1.97(L00) 


1.43 


A2163 


0.2030 


13.29+»«^(W00) 


23.435 ± 1.50(B04) 


82.021 ±5.24 


18.44 ±0.24(FF01) 


1.86 


A2744 


0.3080 


8.65t°-g(Z04) 


13.061 ± 2.44(B04) 


37.315 ±6.97 


17.16 ±1.71(GFG01) 


1.64 


A2219 


0.2280 


9.52t0-55(W00) 


12.732 ±0.98(E98) 


40.293 ±4.34 


12.23 ±0.59(B03) 


1.56 


CL0016-I-16 


0.5545 


9.131°-^^ (WOO) 


18.829 ± 1.88(T96) 


51.626 ±5.16 


6.74 ± 0.67(GF00) 


0.79 


A1914 


0.1712 


10.53i;^'^J(W00) 


10.710 ±1.02(E96) 


33.738 ±3.21 


5.21±0.24(B03) 


1.18 


A665 


0.1816 


8.40+*-g(M96) 


9.836 ±0.98(E98) 


25.130 ±3.92 


3.98 ± 0.39(GF00) 


1.69 


A520 


0.2010 


7.84t°;g(m) 


8.830 ± 0.79(E98) 


22.841 ±5.14 


3.91±0.39(GFG01) 


1.00 


A2254 


0.1780 


7.50«:g(e) 


4.319 ±0.26(E96) 


11.076 ±0.66 


2.94±0.29(GFG01) 


0.86 


A2256 


0.0581 


6.90ton(WOO) 


3.814 ±0.16(E96) 


9.535 ±0.42 


0.24±0.02(F00) 


0.85 


A773 


0.2170 


8.39+0'g(m) 


8.097 ± 0.65(E98) 


21.728 ±3.62 


1.73±0.17(GFG01) 


1.14 


A545 


0.1530 


5.50+f'^0(D93) 


5.732 ±0.50(B04) 


12.608 ±1.10 


1.48 ± 0.06(B03) 


0.82 


A2319 


0.0559 


8-84l:°'^9(M98) 


7.403 ±0.41(E96) 


20.730 ±1.14 


1.12±0.11(F00) 


1.01 


A1300 


0.3071 


9.42«:f5(ml) 


14.114 ±2.08(B04) 


33.870 ± 4.98 


6.09±0.61(F00) 


0.86 


A1656 


0.0231 


8.2lt°'l^(H93) 


3.772 ±0.10(E96) 


10.182 ±0.26 


0.72+0-0^ (m2) 


0.78 


A2255 


0.0808 


6.87to'2o(WOO) 


2.646 ±0.12(E96) 


6.611 ±0.30 


0.89 ± 0.05(G04) 


0.88 


A754 


0.0535 


9.38^°'^:^ (WOO) 


4.314 ±0.33(E96) 


12.946 ±0.98 


1.08±0.06(B03) 


0.96 



Table 2. Parameters of the /3-fit and cluster mass estimated for the 16 galaxy clusters with GRHs for which /3-fits are avaiable. Col.(l): 
Cluster name. Col. (2): /3-parameter value with Icr error. Col. (3): Core radius in units of hZ„ kpc and corresponding uncertainty. Col. (4): 
Virial mass and is uncertainty in units of /i^q 10*^ ^Q- Col. (5): Virial radius in units of /i^q kpc. Col. (6): Mass estimated inside the 
core radius in units of h^^ 10^^ Mq. Ref. for the (data) source in brackets: (a) Markevitch et. al 2002 (Chandra); (b) RB02 (ROSAT for 
/3-fit and T as in table 1); (c) Govoni ct al. 2001 (ROSAT); (d) Ettori & Fabian 1999 (ROSAT); (c) Ettori ct. al 2004 (Chandra); (f) 
Feretti 2004 (Einstein); (g) Lcmonon ct al. 1997 (ROSAT). 
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13 


Tc 


M„ 


Rv 


Mc 


name 




[kpc /.-;] 


[10*5 Mq] 


[kpc h-^] 


[1013 Mq] 


lE50657-558(a) 


0.70 ±0.07 


179 ± 18 


3.43 ± 0.38 


3301 


9.50 ±1.40 


A2163 (b) 


0.80 ± 0.03 


371 ±21 


4.32 ±0.26 


3766 


22.00 ±1.84 


A2744 (c) 


1.00 ± 0.08 


458 ± 46 


2.87 ±0.26 


3096 


22.10 ±2.96 


A2219 (d) 


0.79 ± 0.08 


343 ± 34 


2.52 ±0.28 


3104 


14.40 ±2.16 


CL0016+16 (e) 


0.68 ±0.01 


237 ± 80 


1.47 ±0.05 


2166 


8.27 ±0.38 


A1914 (b) 


0.75 ±0.02 


165 ± 80 


2.90 ±0.15 


3356 


7.28 ±0.51 


A665 (f) 


0.74 ±0.07 


350 ± 35 


1.97 ±0.30 


2933 


12.10 ±2.20 


A520 (c) 


0.87 ±0.08 


382 ± 50 


2.22 ±0.25 


3018 


14.50 ±2.51 


A2256 (b) 


0.91 ±0.05 


419 ± 28 


2.23 ±0.13 


3281 


14.70 ±1.28 


A773 (c) 


0.63 ± 0.07 


160 ± 27 


1.52 ±0.19 


2636 


4.72 ±0.98 


A545 (d) 


0.82 ± 0.08 


286 ± 29 


1.25 ±0.84 


2562 


7.20 ±4.89 


A2319 (b) 


0.59 ±0.01 


204 ± 10 


1.71 ±0.07 


3009 


5.95 ±0.38 


A1300 (g) 


0.64 ±0.01 


171 ± 80 


1.71 ±0.06 


2609 


5.76 ±0.33 


A1656 (b) 


0.65 ±0.02 


246 ± 15 


1.83 ±0.07 


3136 


7.38 ±0.53 


A2255 (b) 


0.80 ±0.05 


419 ± 28 


1.76 ±0.12 


2996 


12.80 ±1.22 


A754 (b) 


0.70 ± 0.03 


171 ± 12 


2.42 ±0.11 


3379 


6.25 ±0.52 
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Figure 3. Panel a): correlation between tlie radio power at 1.4 GHz and the temperature for the GRHs; Panel b): correlation between 
the radio power at 1.4 GHz and the X-ray temperature for a total sample of 24 cluster with a GRHs or with a smaller size (^ 200 — 700 
kpc h-^). 



which is just visible in the NVSS, versus the X-ray lumi- 
nosity of the hosting clusters. We note that all GRHs have 
a radio surface brightness which is well above that of A3562 
(with the exception of A2256 not visible in the NVSS). The 
fact that clusters in the redshift range z ^ 0.15 — 0.3 have 
similar radio brightnesses indicates that the correlation in 
FigQis not driven by the radio surface limit. Indeed, these 
clusters have Lx ^ 3 • 10** — 3 • 10*^ erg/s and range over 
more than one order of magnitude in radio power, whereas 
the effect of brightness dimming, due to the small z-range, is 
limited to within a factor ~ 1.6. We also note the presence 
of a trend, with the average radio brightness increasing with 
X-ray luminosity (see also Feretti 2004), which further sup- 
ports the notion that the correlation in Fig^at high lumi- 
nosities {Lx ^ 10"*^ erg/s) is not driven by selection effects. 
Furthermore, relatively deep upper limits for non radio-halo 
clusters are now available and in some cases lie well below 
(a factor of > 10) the trend in Fig0 at X-ray luminosities 
> 5 ■ 10** erg/s (Dolag 2006). 

On the other hand one may argue that the NVSS tends 
to select only the most powerful GRHs associated with 
Lx ^ 5 • 10** erg/s clusters. To evaluate the effect of a pos- 
sible bias at these luminosities, we perform a fit by consid- 
ering only GRH clusters with Lx ^ 5 ■ 10** and find a slope 
2.22 ± 0.36 which is consistent within 1 a with Eq. Q 

Thus we conclude that despite the poor statistics, the 
derived correlations (Eq0 stands on sound observational 
basis. 



2.2 Radio Power— ICM temperature correlation 

We also investigate the correlation between the radio power 
at 1.4 GHz and the X-ray ICM temperature. A Pi. 4 - T 
correlation was first noted by Liang et al.(1999) and Co- 



lafrancesco (1999); with a sample of only 8 radio halos 
the last author obtained a steep trend of the form Pi. 4 ex 
T ' -2.08, In Fig. 1^ we report the best fit for our sample. 
The fit has been performed using the form: 



log 



Pi.^ 



3.16 -1024/1: 



Af + bf log 



T 



8keV 



(2) 



and best fit parameters are: Af = —0.390 ± 0.139 and 
bf = 9.83 ± 4.92. We note that the observed P1.4 - T 
correlation is very steep, it seems rather a "wall" than a 
correlation and it is dominated by the large errors of the 
cluster temperatures avaiable to date. In order to test the 
strength of this correlation we try to increase the sample 
by including also 7 additional clusters with smaller (size 
~ 200 — 700 kpc h'^g) radio halos. In Fig. I^J we report the 
best-fit Pi. 4 — T obtained for the extended sample, which 
has a slope bf — 6.40 ± 1.64. Given the large uncertainties 
we note that the two correlations are consistent at the la 
level and, in addition, the value of the lower allowed bound 
of the two slopes is almost the same. 



2.3 Radio Power — virial mass correlation 

The most important correlation for our study is that be- 
tween the virial mass (M„) of a cluster and the radio power 
at 1.4 GHz. This correlation is indeed extensively used in the 
calculations of the RHLFs and number counts (Sec. 7) and in 
constraining the values of the magnetic field in galaxy clus- 
ters to be used in our calculations (Sec. 3). On the other 
hand, this is also the most difficult correlation to derive 
since it is very difficult to measure the cluster mass. Gov- 
oni et al. (2001) first obtained a correlation between the 
radio power and cluster gravitational mass (within 3 h'^g 



6 R. Cassano, G. Brunetti, G. Setti 



46 - 



45 - 



"44 - 



b,= l,47±0.08 
b,= 1.41±0.10 



43 




HIFLUGCS _ 
Radio Halos_ 



42 



13.5 14 14.5 15 

Log(MyMj 



15.5 



Figure 4. Correlation between the X-ray luminosity [0.1-2.4] keV 
and the virial cluster mass: for the HIFLUGCS sample (black 
points) plus the 16 clusters with GRHs (red points, excluding 
A2254 for which no information on the /3-model are available) 
(solid line) and for the HIFLUGCS sample alone (dashed line). 



Mpc radius) estimated from the surface brightness profile of 
the X-ray image using 6 radio halo clusters. This correlation 
was confirmed by Feretti (2003) who extended the sample 
to 10 cluster radio halos and obtained a best fit of the form 
Pi. 4 oc M^'^ , where M is, again, the gravitational mass com- 
puted within 3 h'^g Mpc from the cluster center. However 
it should be pointed out that while the X-ray mass deter- 
mination method gives good results in relaxed clusters, it 
may fail in the case of merging clusters (e.g., Evrard et al. 
1996, Roettiger et al. 1996; Schindler 1996). This is because 
the merger may cause substantial deviation from hydrostic 
equilibrium and spherical symmetry. As a result, the masses 
in merging clusters can be either overestimated (up to twice 
the true mass in the presence of shocks) or underestimated 
(since substructures tend to flatten the average density pro- 
file giving an underestimation of the order of 50 % of the 
true mass; see Schindler 2002). In addition, if the temper- 
ature systematically decreases with increasing radius, then 
the isothermal assumption leads to an overestimation of the 
cluster mass of about 30% at about six core radii (Marke- 
vitchet al. 1998). 

The effect of the scattering produced by all these un- 
certainties can hopefully be reduced by making use of large 
cluster samples. Thus, we choose to obtain the Piaghz — M^ 
correlation by combining the L^ — M„ correlation, obtained 
for a large statistical sample of galaxy clusters, with the 
Pi A — Lx correlation previously derived (Eq0 Fig. 0. We 
use a complete sample of the X-ray-brightest clusters (HI- 
FLUGCS, the Highest X-ray FLUx Galaxy Cluster Sample) 
compiled by Reiprich & Bohringer (2002) (hereafter RB02). 
We use this sample of luminous clusters {L^ ~ 10^^ — 10''^erg 
s~^) since it is large and homogeneously studied. It consists 
of 63 bright clusters with galactic latitude \hii\ > 20°, flux 



/x(0.1 - 2AkeV) > 2 X 10"" ergs s"^ cm"^ and it covers 
about 2/3 of the whole sky. 

The clusters have been reanalyzed in detail by RB02 using 
mainly ROSAT PSPC pointed observations. RB02 report 
the value of /3 and core radius, r-c, for all the 63 clusters 
obtained by fitting the surface brightness profile of the X-ray 
image with a standard /3-model. Then under the assumption 
that the intracluster gas is in hydrostatic equilibrium and 
isothermal (using the ideal gas equations), the gravitational 
cluster mass within a radius r is given by (e.g., Sarazin 1986): 



Mtot{< r) = 



ZKbTt-^P 



Vr2+rij 



(3) 



finipG \r^ + ri 

Eql^gives the total (dark matter plus gas) mass of the clus- 
ter as a function of radius; then one must define a physically 
meaningful radius to compare masses of different clusters. 
It is frequently convenient to use r2oo meaning the radius 
within which the mean total mass density is 200 times the 
critical density of the universe at the cluster redshift. This 
is because M200, the mass contained within 7-200, is usually 
taken as a good approximation of the virial mass since in 
the spherical collapse model the ratio between the average 
density within the virial radius and the mean cosmic den- 
sity at redshift z is Ac = 187r'^ ~ 178 independent of the 
redshift for Qm = 1 (e.g., Lacey & Cole 1993). In general, 
the value of Ac depends on the adopted cosmology. In the 
ACDM cosmology Ac is given by (Kitayama & Suto 1996): 

n0.9052n 



Ac(2) = 187r^(l + 0.4093(.^(2)"-'"""), 



where u)(yz) = ^j{z) 
f^m.ofl 4 



1 with: 



(4) 



(5) 



"■'■^^^^ an,0(l+z)3 + f^A' 

Thus, by using Eq|21we calculate the virial radius, R^, as 
the radius at which the ratio between the average density in 
the cluster and the mean cosmic density at the redshift of 
the cluster is given by Ac(z) (Eq^J- The virial mass, M„, 
and the virial radius are thus related by: 

-|l/3 

(6) 



Rv 



3M„ 



47rAc(z)pm(«) 



where p„(z) = 2.78 x 10" VLm,o (1 + ^)^ h^ M^Mpc'^ is the 
mean mass density of the universe at redshift z. 

We estimate the virial mass in the ACDM cosmology 
for the 63 clusters of the HIFLUGCS sample using Eq. |3] 
the fit parameters (/3 and Tc corrected for a ACDM cosmol- 
ogy) and the temperature T are given in RB02. We have 
searched in the literature for /3-fit parameters and T of the 
clusters with GRHs (ref. in Tab. I^J in order to estimate 
M„ also for these clusters. Since some clusters of the HI- 
FLUGCS sample are also in our sample, we note that in 
the majority of these cases the fits to the mass profile (and 
T) given in RB02 leads to a virial mass which is consis- 
tent at Icr level with the mass derived by making use of the 
parameters obtained from more recent observations in the 
literature (given in Tabs. 1,2). The L^ — M„ distribution of 
the combined sample is reported in Fig. |1| The presence of 
a relatively large dispersion indicates the difficulty in esti- 
mating the virial masses of the single objects and confirms 
the need of large samples in these studies. We note that the 
statistical distribution of clusters with GRHs is not different 
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from that of the HIFLUGCS sample. On the other hand, we 
note that clusters with known GRHs span a range in mass 
comparable to the mass-dispersion in the HIFLUGCS sam- 
ple which is due to the different dynamical status of clusters 
in the sample and to the uncertainties in the measurements. 
This further strengthens the need of the approach followed 
in this Section, since a Lx (or Pi.4)-M„ fit based on GRHs 
alone would be affected by large uncertainties. 

In order to better sample the region of higher X-ray 
luminosities and masses (typical of clusters with GRHs) , we 
compute the L^-M^ fit by combining the HIFLUGCS with 
the radio-halo sample. The fit has been performed using the 
form: 



log 



f 0*4 h. 



-1 ergs 



Af+bf log ( 



M„ 



3.16 X 1014 /ifoM© 



(7) 



The best fit values of the parameters are: Af — —0.229 ± 
0.051 and bf = 1.47±0.08 (6/ = 1.41±0.10 is obtained with 
HIFLUGCS sample only). 

In order to derive the Piaghz — M^ correlation for 
GRHs, we combine EqsQandQand find : 



4 ^ 



3 ^ 



2 - 




log 



Pi.. 



3.16- 1024 /ii 



'70 Hz 

(0.814 ±0.147) 



= (2.9 ± 0.4) log 



M„ 



1015 h-Q Mq 



(8) 



Our PiAGHz — Mv correlation is slightly steeper than 
that obtained with 10 clusters by Feretti (2003) {Piaghz oc 
M^'^), which, however, was derived in an EdS cosmology 
by considering the mass within 3 /i^q Mpc from the cluster 
centers, and not the virial mass. 



Figure 5. Expected slope of the Pi. 4 — My correlation as a 
function of the magnetic field intensity in a cluster of mass 
< M >= 1.6 X IQI^Mq. The calculations are obtained for 
b=0.5,0.6,0.7,0.8,0.9,l,1.2,1.3,1.5 and 1.7 (from bottom to top); 
Ml = 1.1 X 10^5 Mq and M2 = 2.5 X lO^^ Mq are adopted. 
The continuous lines are for F ~ 0.67 and the dashed lines are 
for r ~ 0.56. The two horizontal lines mark the 1 a value of the 
observed slope. 



3 EXPECTED CORRELATIONS AND 
MAGNETIC FIELD CONSTRAINTS 

The main goal of this Section is to extract the values of the 
physical parameters to be used in the model calculations of 
Sec. 4-6. The region of the physical parameters (in particular 
of B) is constrained by comparing the model expected and 
observed trends of the synchrotron power of GRHs with the 
mass (and temperature) of the parent clusters. As already 
discussed (Sec. 2) it is unlikely that the observed correla- 
tions are driven by selection effects; however, it cannot be 
excluded that the detailed shape and scatter of these cor- 
relations might somewhat change with improved statistics, 
especially at low X-ray luminosities. 



3.1 Radio power— cluster mass correlation 

Cassano & Brunetti (2005) derived an expected trend be- 
tween the bolometric radio power, Pr, and the virial clus- 
ter's mass and/or temperature. In the case of the GRHs, 
the mergers which mainly contribute to the injection of tur- 
bulence in the ICM are those with Vs ^ Rh, r^ being the 
stripping radius of the infalling sub-cluster (see Sec. 6 in 
CB05). It can be shown that, as a first approximation, the 
expected scaling Pr — M„ is given by: 




PflOC 



(S' + BLJ 



[mG] 



Figure 6. Expected slope of the Pi. 4 — T correlation as a 
function of the magnetic field intensity in a cluster with tem- 
perature < T >= 8 keV. The calculations are obtained for 
b=0.5,0.6,0.7,0.8,0.9,l,1.2,1.3,1.5 and 1.7 (from bottom to top); 
Ti = 6 keV and T2 = 10 keV are adopted. The continuous lines 
are for V ~ 0.67 and the dashed lines are for F ~ 0.56. The two 
horizontal lines mark the 1 a value of the observed slope. 



(9) 
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where B is the rms magnetic field strength in the radio 
halo volume (particle pitch angle isotropitazion is assumed), 
Bcmb = 3.2(1 + z) ^G is the equivalent magnetic field 
strength of the CMB and n^ is the number density of rela- 
tivistic electrons in the volume of the GRH. The parameter 
r is defined by T oc M^ \ we consider F ~ 2/3 (virial scaling) 
and r ~ 0.56 (e.g.; Nevalainen et al. 2000). 
In this paper we release the assumption adopted in CB05 
of a magnetic field independent of cluster mass and assume 
that the rms field in the emitting volume scales as B = 
B<M>{M/ < M >)'', with 6 > and B<m> the value of 
the rms magnetic field associated to a cluster with mass 
equal to the mean mass < M > of the clusters sample. A 
scaling of the magnetic field intensity with the cluster mass 
is indeed found in numerical cosmological MHD simulations 
(e.g. Dolag et al. 2002, 2004). Dolag et al. (2002) found a 
scaling B cc T^ that would mean B oc A/^'^^ assuming the 
virial scaling or B oc M^'^^ for F ~ 0.56. 
We assume that the number density of the relativistic elec- 
trons in galaxy clusters, Tie, does not depend on cluster mass. 
This is because there is no straightforward physical reason to 
believe that this value should scale systematically with A/„, 
and since only a relatively fast scaling of n^ with mass would 
significantly affect the radio power - mass trend (Eq. |5J . It 
is indeed more likely that n^ may change from cluster to 
cluster, but that the major effect would simply be to drive 
some scattering on the Pr — Mv trend (Eq.lSJ. 
Given these assumptions Eq. [^becomes: 



PflOC 



Mfr^ Bt 



(M„/ < M >f 



{B^ 



<M> 



{M,/<M>Y^ + Bl^,Y 



(10) 



which has two asymptotic behaviors: Pr oc M„ + for 
-B<Af> << Bcmb and Pr oc M^"'""^*' for B<Af> >> Bcmb- 
The observed correlations derived in Sect. 2 involve the 
monochromatic radio power at 1.4 GHz. How this monochro- 
matic radio power can be scaled to Pr depends on the spec- 
trum of radio halos. In the context of particle acceleration 
models (e.g., Brunetti et al. 2001, Ohno et al 2002, Kuo et 
al. 2003) the spectrum of radio halos is given by the su- 
perposition of spectra emitted from regions in the emitting 
volume with different magnetic field strenghts. It is expected 
to reach a peack at vi, and then gradually drop as a power- 
law which should further steepen at higher frequencies. The 
break frequency can be expressed as a function of the cluster 
mass and of the rms field B in the emitting volume(CB05): 



Vi oc M 



Brit 



iB' + B!^X 



(11) 



If we adopt a power-law spectrum extending from the fre- 
quency of the peak to a few GHz, P(i') oc z^"", Pr and the 
monochromatic radio power at a fixed frequency i/o iyo ^ Vb) 
scale as P{vo)/Pr oc ( — )°~'^. This depends on the cluster 
mass (Ea lllH : 



Pi^o) 



fB3,,-.(M„/<M>)2'' + B2 



,2{a-l) 



(12) 



^<M>V 



cmb 1 



thus in the case B « Bcmb one has P{vo)/Pr oc 

M \{a-l)(2-T-3b) 



\<M> I 

P{vo)IPr « (<^ 
B « Bcmb the P{i^o 



^cmb 

which means that for 
-M trend is steeper than the Pr — M, 
while the opposite happens in the case B » Bcmb (the two 



Table 3. Values of am and r]t derived for relevant sets of b, 
B^mylfiG] parameters. 



b 


B<Af > [mC] 


«A/ 


Vinin 


f)max 


1.7 


3.0 


2.5 


0.19 


0.2 


1.7 


2.2 


3.22 


0.17 


0.2 


1.5 


1.9 


3.3 


0.15 


0.2 


1.3 


2.25 


2.84 


0.15 


0.2 


1.0 


1.55 


2.96 


0.16 


0.21 


1.0 


0.45 


3.3 


0.29 


0.33 


0.9 


0.18 


3.23 


0.39 


0.44 


0.6 


0.2 


2.63 


0.38 


0.44 



scaling should be equal for continuity for B ~ Bcmb)- On the 
other hand, the trends of P{vo)/ Pr with the cluster mass in 
massive galaxy clusters is rather weak because the observed 
radio spectral index between 327-1400 MHz is a ~ 1.2 (e.g., 
Feretti 2003) and because B in the most massive objects 
is probably close to Bcmb (Sec. 3. 3, Fig|7| Govoni & Feretti 
2004). Thus, in order to compare the model expectations 
with the observations, we will safely assume the same scal- 
ing for monochromatic and total radio power. 

In order to have a prompt comparison with observa- 
tions we calculate the slope um of the Pi, 4 — Ad correlation 
between two points as: 



CtM 



log(Pl/P2) 

log(Afi/M2) 



(13) 



Ea ll3l can be compared with the observed slope to constrain 
the value of the magnetic field and of the slope, 6, of the 
scaling between B and the cluster mass. The M\ and A/2 
values give the representative mass range spanned by the 
bulk of clusters with GRHs, while Bcmb should be calculated 
at the mean redshift of our sample (< z >~ 0.19). We point 
out that given B^m> and b, the values of B are fixed for all 
the values of the masses of the clusters in our sample. 

In Fig|S] we report the expected slope om (Eq. I13II as 
a function of B^m>- The different curves are obtained for 
different scaling-laws of the magnetic field with the cluster 
mass (6 — 0.5 to 1.7, see caption). Dashed lines refer to 
F ~ 0.56 and solid lines to the virial case. The two blue 
horizontal lines (Fig|SJ indicate the range of the observed 
slope {au = 2.9 ±0.4). 

Fig.|S]shows that there are values of B^m> and b for which 
the expected slope is consistent with the observed one. As 
a first result we find that with increasing b the values of 
P<Af> should increase in order to match the observations 
(for example, b ~ 0.6 requires B<^m> ~ 0.2 — 1.4 /iG 
while b ^ 1.7 requires B<^m> ~ 2 — 3 /^G). Finally, the 
asymptotic behavior of Eg 1101 combined with the observed 
correlation (Eq. |HJ allows to immediately constrain b: for 
B<A/> << Bcmb one has 0.58(0.53) < b < 0.98(0.93) for 
the virial (non-virial) case, whereas in the case of B<^m> » 
Bcmb the model expectations cannot be reconciled with the 
observations. 
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Figure 7. The region in the plane (B<j\f>,b) allowed from the 
observed Pi. 4 — My and Pi, 4 — T correlations is reported as a 
shadowed area; < M >= 1.6 X lO^^Af0. The dashed line indi- 
cate the upper bound of the allowed region obtained considering 
only the Pi. 4 — M„ correlation. The coloured points indicate the 
relevant configurations of the parameters used in the statistical 
calculations in Sec. 4-6 (Tab.UJ. The vertical arrows indicate the 
IC limits on B. 



3.2 Radio power— cluster temperature correlation 

Since the temperature is related to the cluster mass, the 
radio power - mass correlation also implies a correlation 
between synchrotron radio power and cluster temperature. 
Thus, in order to maximize the observational constraints, an 
analysis similar to that of Sect. 3.1 can also be done for the 
radio power - temperature correlation {Pr ~T). Combining 
Eq.[Tniwith the M - T scaling law (T ot M^/^ for the virial 
case and T ex M°-^^) one has: 



PflOC 



TT 



-1 d2 



(T/ < r >) 



{B'' 



<M> 



(T/ < r>)2bT 



(14) 



where br ~ b/T with F ~ 2/3 (virial case) or F ~ 0.56 (non- 
virial case). The asymptotic behaviors of Ea. ll4l are given by 



« Be 



and Pr oc M„ 



Pr oc T^/r-i+^fcT ^B^M> 

{B<M> >> Bcmb)- 

As in Sec. 3.1, here we can adopt the same scaling with T 
for both Pr and Pi. 4 and compare the values of the expected 
slope with those of the observed one. We can calculate the 
slope ot of the Pi. 4 — T correlation between two points as: 



CtT 



iog(Pi/P2: 

iog(Ti/r2: 



(15) 



where Ti and T2 define the interval of temperature of our 
sample, < T >= 8 keV is the mean temperature, and Bcmb 
is evaluated at < z >~ 0.19. In Fig.|S|we report the slope ar 
of the P1.4 — T correlation as a function of the magnetic field 
strength in a average cluster, S<m> • The different curves are 
obtained for different scaling-laws of the cluster magnetic 
fields with mass (i.e., temperature) (b=0.5 to 1.7). Dashed 



Figure 8. ProbabiUty to form GRHs at 0.05 ^ z ^ 0.15 in the 
observed mass bin I: 0.95-1.9 X lO^^M© and at 0.05 ^ ^ ^ 0.2 in 
bin II: 1.9 — 3.8 X 10^^ Mq as a function oirjt- The calculations are 
reported for the following representative cases: b = 1.7, B^my = 
3.0/jG (blue points); b = 1.0, B^m> = 1.55^G (black points); 
b = 0.9, B<A/> = 0.18/^G (cyan points) and b = 0.6, P<j\/> = 
0.2/iG (green points). The bottom shadowed region marks the 
observed probability for GRHs in the mass bin I while the top 
shadowed region marks that in the mass bin II. The values of the 
observed probabilities are obtained by combining the results from 
Giovannini et al. 1999, Giovannini & Feretti 2000, and Feretti 
2002. The observed probabilities for the bin I are calculated up 
to 2 ^ 0.15 to minimize the effect due to the incompleteness of 
the X— ray and radio catalogs used by these authors. 



lines are for F ~ 0.65 and continuous lines are for the virial 
case. 

The horizontal blue lines mark the lower limit qt — 
4.76 and the upper limit ar — 8.05 of the observed cor- 
relation. Fig. |S| shows that there is a range of values of the 
parameters (P<i\/>, 6) for which the model is consistent with 
the observed slope. The relevant point is that, similarly to 
the case of the P1.4 — M correlations, also in this case values 
of B<A/> >> Bcmb cannot be reconciled with observations: 
a clear upper boundary at P < 3^G is obtained for B^m>- 



3.3 Constraining the magnetic field 

We combine the results obtained from the observed correla- 
tions (both Pi, 4 — Ml, and Pi. 4 — T) and the model expected 
trends to selects the allowed region of the (B<a/>, b) pa- 
rameters. We consider the slope of the P1.4 — T correlation 
ut — 6.4 ± 1.64 as derived for the extended sample of 24 
galaxy clusters with giant and small radio halos. This is be- 
cause what is important here is the allowed lower bound 
of the values of ut which does not depend on the adopted 
sample (Sect. 2.2). 

In FigEJwe report the region of the plane (P<m>, b) allowed 
by the observed slopes at la level. The lower bound of the 
(P<M>,b) region is due to the Pi. 4 — M„ correlation while 
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Figure 9. a) Occurrence of GRHs as a function of the cluster mass in three redshift bins:0-0.1 (black line), 0.2-0. 3 (blue line), 0.4-0.5 
(green line), b) Occurrence of GRHs as a function of redshift in two mass bins: [1-2]x101^^Mq (cyan line) and [2-4.5] xlO'^^M© (blue 
line). The calculation have been performed assuming: b=1.7, iJ<jv/> = S.O^G, rjt = 0.2 in both panels. 



the upper bound is mostly due to the Pi. 4 — T correlation 
which is poorly constrained because of the very large sta- 
tistical errors. This bound is however also limited by the 
Pi. 4 — Mv correlation (FigQ dashed line). 

An additional limit on B^m>, also reported in Fig|7] 
(vertical arrows), can be obtained from inverse Compton 
(IC) arguments. Indeed a lower bound to the magnetic field 
strength can be inferred in order to not overproduce, via IC 
scattering of the photons of the CMB radiation, the hard-X 
ray excess fluxes observed up to now in a few clusters (e.g., 
Rephaeh & Gruber 2003, Fusco-Femiano et al 2003). In this 
case the value of the mean magnetic field intensity in the 
cluster volume can be estimated from the ratio between the 
hard-X ray and radio emission. The resulting value of the 
magnetic field should be considered as a lower limit because 
the IC emission may come from more external region with 
respect to the synchrotron emission (e.g., Brunetti et al. 
2001, Kuo et al. 2003, Colafrancesco et al. 2005b) and also 
because, in principle, additional mechanisms may contribute 
to the hard-X ray fluxes (e.g., Fusco-Femiano et al. 2003). 
One of the best studied cases is that of the Coma cluster 
for which an average magnetic field intensity of the order 
of Bic — 0.2 /iG was derived (Fusco-Femiano et al. 2004). 
As a first approximation we can use this value to obtain 
the lower bound of B for each cluster mass from the scaling 
B = P<Af>(M/<A/>)^ 

The resulting (B<Af>,b) region spans a wide range of 
values of B and b. An inspection of Fig[7| immediately iden- 
tifies two allowed regimes: a super-linear scaling {b > 1) with 
relatively high values of B and a sub-linear scaling (6 < 1) 
with lower values of B. 

All the calculations we will report in the following sec- 
tions are carried out by assuming representative values of 



(P<A/>,6) inside the constrained region (Fig. |7| coloured 
filled dots and TabOJ. 



4 PROBABILITY TO FORM GIANT RADIO 
HALOS 

4.1 Probability of radio halos and constraining rjt 

In this Section we derive the probability with cluster mass to 
find GRHs in the redshift range 2=0-0.2. The byproduct of 
the Section is to calibrate the model by requiring that the 
expected fraction of cluster with GRHs is consistent with 
the observational constraints. This allows to select a range 
of values of the parameter r/t, which is the ratio between 
the energy injected in the form of magnetosonic waves and 
the PdV work done by the infalling subclusters in passing 
through the most massive one. i]t is a free parameter in our 
calculations since the fraction of the energy which goes into 
the form of compressible modes is likely to depend on the 
details of the driving force. 

In the conservative case of solenoidal forcing (and beta of 
plasma >> 1) this fraction is expected to scale with M^ TZe 
(with A^s < 1, the turbulent Mach number) for Ml Tie. < 
10 and with a flatter slope for larger values (Bertoglio et al. 
2001). Assuming a Reynolds number (at the injection scale. 
I.e., hundreds of Kpc) in hot and magnetized galaxy clusters 
Tie X^IO^ (see discussion in Lazarian 2006; Brunetti 2006) 
and a turbulent energy of the order of ~ 20% of the thermal 
energy (CB05), from Fig. 8 in Bertoglio et al. (2001) one 
finds a reference value rjt ~ 0.1 which may be even larger in 
the case of compressible driving. 

Radio halos are identified with those objects in a syn- 
thetic cluster population with a synchrotron break frequency 
(Eqinj Vb > 200 MHz in a region of 1 Mpc h'^r^ size. In 
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Figure 10. a) Occurrence of GRHs as a function of the cluster mass in three redshift bins: 0-0.1 (black line), 0.2-0. 3 (blue line), 0.4-0.5 
(green line), b) Occurrence of GRHs as a function of redshift in two mass bins: [1-2]x101^^Mq (cyan line) and [2-4.5] xlO'^^M© (blue 
line). The calculation have been performed assuming: b=0.9, i?<A/> = 0.2^G, »?t = 0-42 in both panels. 



CB05 it was assumed that the magnetic field in the radio 
halo volume is independent from the cluster mass and it is 
B ~ 0.5/iG. Then ui, oc M^~^ and consequently mas- 
sive clusters are expected to be favourite in forming GRHs. 
CB05 indeed showed that the expected fraction of clus- 
ters with GRHs naturally shows an abrupt increase with 
cluster mass, and that the observed fractions (20-30 % for 
Af > 2 X lO^'^ Mq clusters, 2-5 % for M ~ 10^^ M© clusters 
and negligible for less massive objects) can be well reconciled 
with the model expectations by assuming rjt ~ 0.24 — 0.34. 

In the present paper we assume that the rms magnetic 
field depends on the cluster mass and this should afi^ect the 
synchrotron break frequency (Eq. IIH and the occurrence of 
GRHs with cluster mass. On the other hand, in Sect. 3.3 we 
have also shown that the comparison between the expected 
and observed trends between radio power and cluster mass 
(and temperature) helps in constraining the range of values 
which can be assigned to the magnetic field in clusters. 
Thus our calculations of the occurrence of GRHs {z ^ 0.2) 
and the selection of the values of rjt necessary to reproduce 
the observations should be performed within the dashed re- 
gion in Fig0 

To calculate the expected probabilities to form radio 
halos we first run a large number. A/", of trees for different 
cluster masses at 2 = 0, ranging from ~ 5 x 10^'^Mq to 
~ 6 X 10^^ Mq. Then we choose different mass bins AM and 
redshift bins Az in which to perform our calculations. Thus, 
for each mass M, we estimate the formation probability of 
GRHs in the mass bin AM and in the redshift bin Az as 
(CB05): 



Z.^,=l *f 



E7=i(ii + *'J 



(16) 



where f„ is the time in the redshift interval Az that the 
cluster spends in the mass bin AM with i/t ^ 200MHz and 
td is the time that the same cluster spends in AM with 
ut < 200 MHz. The total probability of formation of GRHs 
in the mass bin AM and in the redshift bin Az is obtained 
by combining all the contributions (Eq. 1161 weighted with 
the present day mass function of clusters given by the Press 
& Schecther mass function. 

To have a prompt comparison with present observa- 
tional constraints, we calculate the probability to form 
GRHs at 2 < 0.2 in the two observed mass bins: bin I 
([0.95 - 1.9] X lO^^M©) and bin H ([1.9 - 3.8] x IO^^Mq). 

As an example, in Fig.|H|we report these probabilities in 
both bin I and bin H as a function of rjt for three represen- 
tative cases which nicely sample the region in Fig|7| b = 1.7, 
B<M> = 3.0/iG (blue points); h — 1.0, B^m> ~ 1.55^G 
(black points); b — 0.9, B^m> = 0.18/iG (cyan points); 
b = 0.6, B^M^ — 0.2/iG (green points). 

The bottom shadowed region in Fig. |S] marks the ob- 
served probability for GRHs in the mass bin I while the top 
shadowed region marks that in the mass bin H. Fig.|H]shows 
that it is possible to find a range of values of the parame- 
ter r^t for which the theoretical expectations are consistent 
with the observed statistics in both the mass bins. However 
we note that the requirement in terms of energy of the MS 
modes increases with decreasing the magnetic field: it goes 
from rjt '^ 0.15 — 0.2 for intermediate-large values of B up 
to rjt ~ 0.5 at the lower bound of the allowed B strengths. 

The fact that the magnetic field depends on the cluster 
mass is reflected in the different behavior that the various 
selected configurations of parameters may have in Fig. |H|in 
the two mass bins: one configuration of parameters may be 
favoured in a mass bin with respect to another configuration 
but disfavoured in the other mass bin. This is related to the 
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transition from IC dominance [B < Bcmb) to synchrotron 
dominance {B > Bcmb) that occurs in going from the bin 
I to the more massive clusters of bin II. In the case of IC 
dominance an increase of B does not significantly affect the 
particle energy losses, it causes an increase of i^b fEa lllH 
and thus an increase of the probability to have GRHs. On 
the other hand, in the case of synchrotron dominance the 
particle energy losses increase and consequently Ub decreases 
fEa lllll as well as the probability to form GRHs. 

For this reason, given i]t, the ratio between the proba- 
bility to form GRHs in the bin I and in the bin II is expected 
to decrease with increasing b, as larger values of b yield a 
more rapid increase of B with cluster mass (Fig|HJ . 

In Tab 13 we report the maximum and the minimum val- 
ues of rjt {rit,max and rjt,min) for which the model reproduces 
the observed probabilities {la limits) in both the mass bins. 
The results are given for the relevant (_B<a/>,6) configura- 
tions reported in Fig. |7] In agreement with the above dis- 
cussions, one might notice that in the case of IC dominance 
a larger magnetic field implies a smaller energetic request 
(smaller rit,max)- 




IQZ4 10=^== 

P,., [W/Hz h-s] 



10= 



4.2 Probability of radio halos with My 
evolution with z 



and 



In this Section we calculate the expected probability to form 
GRHs with cluster mass without restricting ourselves to the 
mass bins considered by present observations (bin I and bin 
II in Fig. ISJ and calculate the evolution of this probability 
with redshift. In doing these calculations we use the values 
of rjt as constrained in Tab|21within the region (_B<Af>,b) of 
FiglZl S'lid make the viable (and necessary) assumption that 
the value of rjt (i.e., efficiency of turbulence in going into MS 
modes) is constant with redshift. 

A detailed calculation of the acceleration efficiency and 
of the probability to have GRHs requires detailed Monte- 
carlo calculations (see Sec. 6 of CB05) essentially because 
at each redshift the acceleration is driven by MS modes in- 
jected in the ICM from the mergers that the cluster expe- 
rienced in the last few Gyr at that redshift. However, to 
readily understand the model results reported in the follow- 
ing, we may use the simplified formula Eq. (Illll which de- 
scribes the approximate trend of the break frequency with 
cluster mass. The scaling B oc M* adopted in this paper 
implies that the synchrotron losses overcome the IC losses 
first in the more massive objects. Clusters of smaller mass 
in our synthetic populations have B « Bcmb and this im- 
plies (Eqinj n oc M^"'"+* (1 + z)-^ so that the probability 
to form GRHs in these clusters increases with the cluster 
mass (2 — F -f- b > always) and decreases with redshift. In 
the case of more massive clusters the situation may be more 
complicated. Indeed for these clusters there is a value of 
the mass. A/*, for which the cluster magnetic field becomes 
equal to B^mb- For M > M,{z) it is Vb oc M^"""-^' (Eq.ITTJ 
and thus the probability to form GRHs would decrease as 
the mass becomes larger (given the lower bound of the slope 
b as constrained in Fig. Q it is 2 — F — 36 < 0). In these 
cases, at variance with the smaller clusters, the occurrence 
of GRHs with z is only driven by the cosmological evolution 
of the cluster-merger history (which drives the injection of 
turbulence) rather than by the dependence of the IC losses 
with z (at least up to a redshift for which B ~ Bcmb(z)). As 



Figure 11. Expected RHLFs at z ~ 0.05 (coloured lines with 
dots) obtained assuming: b=1.7, B^m> = S.O^G (blue lines: rjt = 
0.2 (solid line) and r]t = 0.19 (dashed line)); b=1.7, B^m> = 
2.2fiG and rjt = 0.2 (magenta line); b=1.5, B^m> = 1.9/iG and 
r]t = 0.2 (red line); b=0.9, B<m> = 0.18/iG and rjt = 0.39 (cyan 
line); b=0.6, -B<ai> = 0.2/iG and rjt = 0.38 (green line); b=1.0, 
S<A/> = 0.45/iG and rjt = 0.33 (black line). For a comparison 
we report the range of Local RHLF obtained by E&R02 (black 
solid thick lines). 



a consequence, the general picture is that going from smaller 
to larger masses, the probability should reach a max value 
around Af, for which B ~ Bcmb{z), and then it should start 
to smoothly decrease. The value of this mass increases with 
z and depends on the scaling law of B with M. It is: 



^^ , ^ ,^ / 3.2(14- z)^\ 



i/b 



(17) 



In order to show in some detail this complex behavior 
in the following we analyze two relevant examples. 



4-2.1 An example with super-linear scaling: large B 

As a first example we focus on the case of a super-linear 
scaling. In Fig. |5] we report the occurrence of GRHs as a 
function of the cluster mass in three redshift bins (panel 
a)) and the occurrence of GRHs as a function of redshift in 
two mass bins (panel b)). These calculations have been per- 
formed using 6 = 1.7 and -B<a/> = 3/xG which are allowed 
from the observed correlations. We adopt rjt ~ 0.2 which is 
in the corresponding range of values obtained in Sec. 5 (see 
Tab. O in order to reproduce the observed probability of 
GRHs at 2: < 0.2. One finds that at lower redshifts (z < 0.1) 
the probability to form GRHs increases with the mass of 
the clusters up to Ai", -- 2 x 10^^ A/q , while for M >M, 
synchrotron losses become dominant and this causes the de- 
crease of the probability for M -^Af, . The mass at which 

B ~ Bcmb{z) increases as (1 -I- z)^''' and this causes the 
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Figure 12. Evolution of RHLFs with redshift. The RHLFs are reported from redshifts 0-0.1 to 0.5-0.6 (curves from top to bottom). 
Calculations are developed for: Panel a) b=1.7, iJ<jv/> = S.O^G, rjt = 0.2, a]\,[ ~ 2.5 and Panel b) b=0.9, i3<i\/> = 0.18 /^G, rjt = 0.39, 
om — 3.23. 



shift with z of the value of the cluster mass at which the 
maximum of the probability is reached. 
FigHt) shows the occurrence of GRHs with z. In the higher 
mass bin (2 ■ 10^^ < M < 4.5 • lO^'') the occurrence in- 
creases up to ~ 0.4 and than start to drop. In this very 
massive clusters the magnetic field is larger than Bcmb{z) at 
any redshift and thus the synchrotron losses are always the 
dominant loss term. The behavior of the probability with 
z in this case is essentially due to the fact that the bulk of 
turbulence in these massive clusters is injected preferentially 
between z ~ 0.2 — 0.5. A different behavior is observed in the 
lower mass bin (10^^ ^ Af ^ 2 • 10^^) where the occurrence 
of GRHs decreases with redshift. This is because clusters 
with these lower masses have always B < Bcmb{z)- 



4-2.2 An example with sub-linear scaling: small B 

As a second example we focus on a sublinear scaling b. In 
Fig. 1101 we report the occurrence of GRHs as a function of 
the cluster mass in three redshift bins (panel a)) and the 
occurrence of GRHs as a function of redshift in two mass 
bins (panel b)). The calculations have been performed using 
b = 0.9 and B^m> ~ 0.2/iG, which are allowed from the 
correlations, and adopting a corresponding rjt — 0.42, which 
is within the range of values obtained in Sec. 5 (see Tab.|^ 
in order to reproduce the observed probability of formation 
of GRHs at redshift z < 0.2. In this case at any redshift 
the probability to form GRHs increases with the mass of 
the clusters. Indeed the magnetic field in these clusters is 
always B « Bcmb{z) (for all redshifts and masses) and the 
IC losses are always the dominant loss term. In addition, as 
expected, in both the considered mass bins the probability 
to form GRHs decreases as a function of redshift, due to the 
increase of the IC losses (Fig. 1101 panel b)). 



5 LUMINOSITY FUNCTIONS OF GIANT 
RADIO HALOS 

In this Section we derive the expected luminosity functions 
of giant radio halos (RHLFs). Calculations for the RHLFs 
are carried out within the (i3<Ai>,b) region of Fig. Q by 
adopting the corresponding values of rjt which allow to 
match the GRH occurrence at z < 0.2. First we use the 
probability P/^^' to form GRHs with the cluster's mass to 
estimate the mass functions of GRHs {dNH{z)/dMdV): 



dNniz) _ dN,i{z) 



X Pf « 



nps X Faz , 



(18) 



dMdV dMdV 
where nps = nps{M,z) is the Press & Schechter (1974) 
mass function whose normalization depends essentially by 
(Tg (present-day rms density fluctuation on a scale of 8h~^ 
Mpc) and Slo; we use as = 0.9 in a fio = 0.3 universe. We 
stress that we use nps since our model is based on the Press 
& Schechter formalism. 
The RHLF is thus given by: 



dNHJz) _ dNajz) I dJ\A_ 
dVdPiA ^ dM dV / dM 



(19) 



dPiA/dM depends on the adopted (B<cm>, b) since each 
allowed configuration in Fig.|2|selects a value of the slope of 
Pi A — Mv (e.g.. Tab. 3) which is consistent (at 1 a) with the 
value of the observed slope obtained with present observa- 
tions (ttM ~ 2.9 ± 0.4; see Sec. 3). In particular from Fig.|^ 
one has that, for a given b, larger values of the magnetic field 
select smaller values of the slope of the Pi a — Mv correlation 
(and viceversa). 

In Fig llll we report the Local RHLFs (number of GRHs 
per comoving Gpc"^ as a function of the radio power) as 
expected from our calculations. The most interesting feature 
in the RHLFs is the presence of a cut-off/fiattening at low 
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Figure 13. Expected RHLFs in 6 redshift bins (as reported in the panels). Calculations are performed by using the following values of 
the parameters : b=1.7, i3<A/> = 3.0/xG (blue lines: r]t = 0.2 (solid lines) and r]t = 0.19 (dashed lines)); b=1.7, B^j^jy = 2.2/iG and 
rjt = 0.2 (magenta lines); b=1.5, B^i^i^ = 1.9/iG and rjt = 0.2 (red lines); b=0.9, B^m> = 0.18^0 and r]t = 0.39 (cyan lines); b=0.6, 
^<u> — 0.2/.tG and rjt = 0.38 (yellow lines); b=1.0, B^m^ = 0A5fiG and rjt = 0.33 (black lines). 
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Figure 14. Number of expected GRHs above a given radio flux 
at 1.4 Ghz from a full sky coverage up to 2 ^ 0.2 (the colour 
code is that of Fig llli . The black points are the data taken from 
Giovannini et al.(1999) and corrected for the incompleteness of 
their sky-coverage (~ 2 7r sr) . 
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Figure 15. Number of expected GRHs from the whole universe 
above a given radio flux at 1.4 GHz. The colour code is the same 
of Figrm 



radio powers. This flattening is a unique feature of particle 
acceleration models since it marks the effect of the decrease 
of the efficiency of the particles acceleration (in 1 Mpc h'^^ 
cube) in the case of the less massive galaxy clusters. We 
stress that this result does not depend on the particular 
choice of the parameters. 



To highlight the result, in Fie llll we also compare our 
RHLFs with the range of Local {RH LF s) EizR (black solid 
lines) reported by Enfilin & Rottgering (2002). These 
(RHLFs) EiiR are obtained by combining the X-ray lumi- 
nosity function of clusters with the radio-X-ray correlation 
for GRHs and assuming that a costant fraction, frh = 1/3, 
of galaxy clusters have GRHs independently from the cluster 
mass (see EnfSlin & Rottgering 2002). 

The most important difference between the two expectations 
is indeed that a low-radio power cut-off does not show up in 
the {RHLFs) EUR in which indeed the bulk of GRHs is ex- 
pected at very low radio powers. The agreement between the 
two Local RHLFs at higher synchrotron powers is essentially 
because the derived occurrence of GRHs in massive objects 
(Sect. 4.2) is in line with the fraction, /^h = 1/3, adopted 
by Enfilin & Rottgering (2002). 

In Fig. 1121 we report the RHLFs expected by our calcu- 
lations in different redshift bins. The calculations are per- 
formed by using two relevant sets of parameters (a super- 
linear and a sub-linear case as given in the caption of 
Fig. I12II allowed from the observed correlations. With in- 
creasing redshift the RHLFs decrease due to the evolution 
of the clusters mass function with z and to the evolution of 
the probability to form GRHs with z. 

Fig. 1121 allows to readily appreciate the different be- 
havior of the RHLFs in the case of a super-linear scaling of 
B with M, & = 1.7, (Fig. EH Panel a)) and of a sub-linear 
scaling, 6 = 0.9 (Fig. 1121 Panel b)): the evolution with red- 
shift in the Panel b) (sub-linear case) is faster than that 
in the Panel a) (super-linear case). This difference is driven 
by the probability to form GRHs as a function of redshift 
in the two cases: in the super-linear case the probability to 
form GRHs does not decrease rapidly with z, while a rapid 
decrease of such a probability is obtained in the sub-linear 
case (see also Figs. |5| 1101 Sect. 6). 

In Fig. ll3l we report the RHLFs obtained by our calcu- 
lations by adopting the selected set of configurations given 
in Tab. |3] (colour code is the same of Fig. 0. The combi- 
nation of these configurations define a bundle of expected 
RHLFs which determines the range of the possible RHLFs. 
All the calculations are performed for the corresponding 
range of values of r]t which allow to be consistent with the 
observed probability to form radio halos at z < 0.2. One 
finds that with increasing redshift the bundle of the RHLFs 
broadens along the nn (P) x P axis. This is again due to the 
different evolutions of the probability to form GRHs with z 
of the super-linear and sub-linear cases. 



6 NUMBER COUNTS OF GIANT RADIO 
HALOS 

In this Section we derive the expected number counts of gi- 
ant radio halos (RHNCs). This will allow us to perform a 
first comparison between the model expectations and the 
counts of GRHs which can be derived from present observa- 
tions, but also to derive expectations for future observations. 
As for the case of the RHLFs, in calculating the RHNCs 
we adopt the configurations of parameters which allow to 
reproduce the observed probabilities of GRHs at z < 0.2. 
However, we point out that the fact that our expectations 
are consistent with the observed probability to form GRHs 
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Figure 16. Expected total number of GRHs above a given radio flux in different redshift bins: panel a) above 5 nijy; panel b) above 30 
mJy. In both panels the colour code is the same of Fig llll 



at «< 0.2 does not imply that they should also be consistent 
with the observed flux distribution of GRHs in the same 
redshift interval. 

Given the RHLFs (diVi/(z)/dPi.4dV) the number of 
GRHs with / > /i,4 is given by: 



Nh{> fi.i) 



dz'{ 



dz" 



J"l.4(/i*4,z') 



dNH{PlA,z' 

dPiA dV 



counts are in excellent agreement with the counts obtained 
from the observations. We note that assuming a superlinear 
scaling of B with cluster mass, up to 30-40 GRHs at z < 
0.2 are expected to be discovered with future deeper radio 
surveys. On the other hand, the number of these GRHs in 
the case of a sublinear scaling should only be a factor of ~ 2 
dPi.4(20) larger than that of presently known halos. 



where dV/dz is the comoving volume element in the ACDM 
cosmology (e.g., Carroll, Press and Turner 1992); the radio 
flux and the radio power are related by Pi. 4 — 47rd|_/i.4 
with di the luminosity distance (where we neglect the K- 
correction since the slope of the spectrum of radio halos is 
close to unity). 

As a first step, we use Eg. l2Ul to calculate the number 
of expected GRHs above a given radio flux at 1.4 Ghz from 
a full sky coverage up to 2:^0.2 and compare the results 
with number counts derived by making use of the present 
day observations ('Fig. 1141 the colour code is that of Fig lllH . 
Calculations in Fig. ll4l are obtained by using the full bundle 
of RHLFs obtained in the previous Section (Fig. 1131 . The 
black points are obtained by making use of the radio data 
from the analysis of the radio survey NVSS by Giovannini 
et al.(1999); normalization of counts is scaled to correct for 
the incompleteness due to the sky-coverage in Giovannini et 
al. (~ 2 7r sr). The NVSS has a la level at 1.4 GHz equal to 
0.45 mjy/beam (beam=45x45 arcsec, Condon et al. 1998). 
By adopting a typical size of GRH of the order of 1 Mpc, the 
surface brightness of the objects which populate the peak of 
the RHLFs (~ 10^'' W/Hz) at z~0.15 is expected to faU 
below the 2a hmit of the NVSS. These GRHs have a flux 
of about 20 mJy, thus below this flux the NVSS becomes 
poorly efficient in catching the bulk of GRHs in the redshift 
bin z=0-0.2 and a fair comparison with observations is not 
possible. For larger fluxes we find that the expected number 



As a second step, we calculate (Fig llSl the whole sky 
number of GRHs expected up 2 = 0.7 (the probability to 
form GRHs at z; > 0.7 is negligible). We note that the 
number counts of GRHs increases down to a radio flux 



of /1.4 ~ 2 — 3 mJy and then flattens due to the strong 
(negative) evolution of the RHLFs fFig. I13II . We note that 
the expected total number of GRHs above 1 mJy at 1.4 
GHz is of the order of ~ 100 depending on the scaling of 
the magnetic fleld with cluster mass. 

Finally we calculate the expected number counts of 
GRHs above a given radio flux in different redshift bins. 
This allows us to catch the redshift at which the bulk of 
GRHs is expected. In Fig. Ilfcil we report the RHNCs inte- 
grated above 5 mJy (Panel a)) and above 30 mJy (Panel 
b)). We note that the bulk of GRHs is expected in the red- 
shift interval 0.1 — 0.3 and this does not strongly depend 
on the flux limit. We note that the relatively high value of 
such redshift range is also due to the presence of the low 
radio power cut-off in the RHLFs which suppresses the ex- 
pected number of low power GRHs. On the other hand, at 
radio fluxes > 30 mJy the contribution from higher redshift 
decreases since the requested radio luminosities at these red- 
shift correspond to masses of the parent clusters which are 
above the high-mass cut-off of the cluster mass function. 
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Figure 17. a) The occurrences of GRHs as a function of the cluster mass in the redshift bins 0-0.1 (solid hues) and 0.4-0.5 (dashed 
lines) are reported for 150 MHz (thick lines) and for 1.4 GHz (thin lines), b) Mass functions of GRHs in the redshift bins 0-0.1 (solid 
lines) and 0.4-0.5 (dashed lines) are reported for 150 MHz (thick lines) and for 1.4 GHz (thin lines), c) Gomparison between the expected 
RHNGs above a given radio flux at 1.4 Ghz (thin lines) and at 150 MHz (thick lines) from a full sky coverage up to z ^ 0.6. 
All the calculations have been performed assuming: b=1.5, i?<M> = 1-9/jG and rjt = 0.2. 



7 TOWARDS LOW RADIO FREQUENCIES: 
MODEL EXPECTATIONS AT 150 MHZ 

Due to their steep radio-spectra, GRHs are ideal targets 
for upcoming low-frequency radio telescopes, such as LO- 
FAR and LWA. In this section we present calculations of 
the statistics of GRHs at 150 MHz derived from the elec- 
tron reacceleration model. 

For simplicity, we present these results only for one set 
of the parameters in the plane (B<a/>,6) (Fig|7J: a super- 
linear case (b=1.5, B<cm> = 1.9/iG) (see Sec. 3). 

First, we calculate the probability to have GRHs at 
^ 150 MHz as a function of the cluster's mass following 
the procedure outlined in Sec. 4.1 and requiring a break 
frequency 2^6^20 MHz to account for the observation fre- 
quency. In Fig ll7b we report the probability to have GRHs 
as a function of virial mass in two redshift bins at 1.4 GHz 
(thin lines) and at 150 MHz (thick lines). As expected, the 
probability at 150 MHz is substantially larger than that cal- 
culated at 1.4 GHz, particularly for higher redshifts and for 
low massive clusters. 

One of the main findings of our work is the presence 
of a cut-off in the RHLFs at low radio powers (see Sec. 5), 
which reflects the drop of the probability to form GRHs as 
the cluster's mass decreases. In Fig |17b we plot the mass 
functions of radio halos (RHMFs) at 1.4 GHz and at 150 
MHz in two redshift bins (see caption of Fig ll7^ . We note 
that the number density of GRHs is increased by only a 
factor -- 2 for Af > 2 • 10^^ Mq, but by more than one 
order of magnitude for M ^ lO'^^ Mq . The most interesting 
feature is again the presence of a low mass cut-off in the 
RHMFs at 150 MHz, which however is shifted by a factor 
~ 2 towards smaller masses with respect to the case at 1.4 
GHz. This is related to the fact that a smaller energy density 
in the form of turbulence is sufficient to boost GRHs at lower 
frequencies, and this allows the formation of GRHs also in 
slightly smaller clusters, which indeed are expected to be 
less turbulent (CB05; see also Vazza et al. 2006). 

Finally, in order to obtain estimates for the RHLFs and 



RHNCs at 150 MHz, we tentatively assume the same Pr — M 
scaling found at 1.4 GHz, scaled at 150 MHz with an average 
spectral index a^ ^ 1.2, and follow the approach outlined 
in Sees. 5 and 6. In Fig llTb we report the expected integral 
number counts of radio halos from a full sky coverage above 
a given radio ffux at 1.4 GHz (thin lines) and at 150 MHz 
(thick lines) up to a redshift z ~ 0.6. The expected number 
of GRHs at 150 MHz are a factor of ~ 10 larger than the 
number expected at 1.4 GHz, with the bulk of GRHs at 
fluxes ^ few mJy. In the near future LOFAR will be able to 
detect diffuse emission on Mpc scale at 150 MHz down to 
these fluxes and this would be sufficient to catch the bulk 
of these GRHs (a more detailed study will be presented in 
a forthcoming paper). 



8 SUMMARY AND DISCUSSION 

The observed correlations between radio and X-ray proper- 
ties of galaxy clusters provide useful tools to constraining 
the physical parameters that are relevant to the reaccelera- 
tion models for the onset of giant radio halos (GRHs). Our 
analysis is based on the calculations of Cassano & Brunetti 
(2005; CB05), which assumes that a seed population of rela- 
tivistic electrons reaccelerated by magnetosonic (MS) waves 
is released in the ICM by relatively recent merger events. To 
this end we have collected from the literature a sample of 17 
GRH clusters for all of which, but one (A2254) , both radio 
and X-ray homogeneous data are available, as summarized 
in Tab.l & 2. Based on the relationships derived in CB05 pa- 
per, we have been able to constrain the (likely) dependence 
of the average magnetic field intensity (B) on the cluster 
mass, under the assumption that B can be parameterized 
as B = B<M>(M/ < M >f (with B<m> the average field 
intensity of a cluster of mean mass < M >= 1.6 x 10^^ Mq 
and b positive). This is an important achievement because 
both the emitted synchrotron spectrum and losses depend 
critically on the field intensity. Following CB05 approach, 
the merger events are obtained in the statistical scenario 
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provided by the extended Press & Schechter formalism that 
describes the hierarchical formation of galaxy clusters. The 
main results of our study can be summarized as follows: 



of the model parameters. 



Probability to form GRHs 



• Observed correlations 

In Sect. 2 we derive the correlations between the radio 
power at 1.4 GHz (Pi. 4) and the X-ray luminosity (0.1- 
2.4 keV), ICM temperature and cluster mass. Most im- 
portant for the purpose of the present investigation is the 
PiA — Mv correlation which has been derived by combin- 
ing the Lx — Mv correlation obtained for a large statistical 
sample of galaxy clusters (the HIFLUGCS sample plus our 
sample) with the P1.4— Lx correlation derived for our sample 
of GRHs. This procedure allows us to avoid the well known 
uncertainties and limits which are introduced in measuring 
the masses of small samples of galaxy clusters, especially in 
the case of merging systems. We find a value of the slope 
om = 2.9 ± 0.4 (P1.4 oc M"^'). A steep correlation of the 
synchrotron luminosity with the ICM temperature is also 
found, although with a large statistical error in the determi- 
nation of the slope : qt = 6.4 ± 1.6 (Pi. 4 tx T"^). 

In Sec. 2 we have also shown that at least in the case 
of high X-ray luminosity clusters {Lx^5 ■ 10^'* erg/s) the 
above trends are unlikely driven by selection effects in the 
present observations. 

• Constraining the magnetic field dependence on the mass 

A correlation between the radio power and the cluster 
virial mass is naturally expected in the framework of electron 
acceleration models. This relationships, discussed in Sec. 3.1 
fEa ll2ll . can reproduce the observed correlation for viable 
values of the physical parameters. For instance, in the case 
B « B„„b, it is P{uo) (X Mi,"<^"'"+'')+'' and the exponent 
agrees with the observed one (om ~ 3) by adopting a typical 
slope of the radio spectrum a = 1 — 1.2 and a sub-linear 
scaling 6 ~ 0.6 — 0.8. 

A systematic comparison of the expected correlations be- 
tween the radio power and the cluster mass with the ob- 
served one (Sect. 3.1 & 2) allows the definition of a permit- 
ted region of the parameters' space (P<A/>,b), where a lower 
bound B<cM> = 0.2 /iG is obtained in order not to overpro- 
duce via the IC scattering of the CMB photons the hard 
X-ray fluxes observed in the direction of a few GRHs (Sect. 
3.3 and Fig. 6). It is found a lower bound at 6 ~ 0.5 — 0.6 and 
that a relatively narrow range of B<m > values is allowed for 
a fixed b. The boundaries of the allowed region, aside from 
the lower bound of B<m>, are essentially sensitive to the 
limits from the Pi. 4 — M„ correlation. 

A super-linear scaling of B with mass, as expected by MHD 
simulations (Dolag et al. 2004) falls within the allowed re- 
gion. 

The values of the average magnetic field intensity in the 
superlinear case are close (slightly smaller) to those ob- 
tained from the Faraday rotation measurements (e.g., Gov- 
oni & Feretti 2004), which, however, generally sample re- 
gions which are even more internally placed than those 
spanned by GRHs. 

Future observations will allow to better constrain the 
radio-X ray correlations and thus to better define the region 



In Sect. 4 we report on extensive calculations aimed at 
constraining rjt, the fraction of the available energy in MS 
waves, which is required to match the observed occurrence 
of GRHs at redshifts z ^ 0.2 (Fig. 7). By adopting a rep- 
resentative sampling of the allowed {B^M>,b) parameter 
space (Fig. 6) we find 0.15 ^ jyt ^ 0.44: the larger values 
are obtained for B^m> approaching the lower bound of the 
allowed region, because of the larger acceleration efficiency 
necessary to boost electrons at higher energies to obtain a 
fixed fraction of clusters with GRHs. 

With an appropriate rjt value for each set of (i3<A/>,b) 
parameters we can calculate the probability of occurrence 
of GRHs at larger redshifts for which observational data 
are not available. This probability depends on the merging 
history of clusters and on the relative importance of the 
synchrotron and IC losses, and shows a somewhat compli- 
cated behavior with cluster mass and redshift. The max- 
imum value of this probability at a given redshift is found 
for a cluster mass M, fEa llTII which mark the transition be- 
tween the Compton and the synchrotron dominated phases. 
In the case of sublinear scaling of the magnetic field with 
cluster mass (b~0.6-0.9) the allowed values of the strength 
of the magnetic field are relatively small (Fig. |7|l, the value 
of M, is large and the IC losses are always dominant 
for the mass range of clusters with known GRHs. As a 
consequence the probability to have GRHs increases with 
cluster mass and decreases with redshift (Fig IIOII . On 
the other hand superlinear scalings (b'-^1.2-1.7) imply 
allowed values of B^m> relatively large (Fig. |7|l, and even 
larger values of the magnetic field for the most massive 
objects. In this case the value Af» falls within the range of 
masses spanned by GRH clusters: the predicted fraction 
of clusters with GRHs increases with mass, then reaches 
a maximum value at about M„ ~ A/*, and finally falls 
down for larger masses (Fig|3Jl. At variance with the case 
of sublinear scaling, in this case the fraction of the most 
massive objects with GRHs is expected to slightly in- 
crease with redshift, at least up to z=0.2-0.4 (Fig|UJ where 
the bulk of turbulence is injected in a ACDM model (CB05) . 

• Luminosity functions (RHLFs) 

In Sect. 5 we report the results of extensive calculations 
following a fair sampling of the (i3<A/>,b) allowed region 
as summarized in Tab. 3; this essentially allows a full cov- 
erage of all possible RHLFs given the present correlations 
at 1(7. We find that, although the large uncertainties in the 
(P<A/>,b) region, the predicted local RHLFs are confined to 
a rather narrow bundle, the most characteristic common fea- 
ture being the presence of a flattening/cut-off at radio pow- 
ers below about 10^* W/Hz at 1.4 GHz (FiglTTJ- The frac- 
tion of GRHs with 1.4 GHz luminosity below ~ 5 x 10^^ W 
Hz~^ h^Q^, a factor of ~ 5 smaller than the luminosity of the 
less powerful GRH (A2256, z=0.0581) known so far, is neg- 
ligible. This characteristic shape of the RHLFs, obtained in 
our paper for the first time, represents a unique prediction 
of particle acceleration models, and does not depend on the 
adopted physical details for the particle acceleration mecha- 
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nism. This is due to the decrease of the efficiency of particle 
acceleration in the case of less massive clusters which is re- 
lated to three major reasons (see CB05): 

i) smaller clusters are less turbulent than larger ones 
since the turbulent energy is expected to scale with the 
thermal one (CB05; see also Vazza et al. 2006); 

ii) turbulence is typically injected in large Mpc regions 
in more massive clusters and thus these are favoured for 
the formation of GRHs (CB05); 

Hi) since in the present paper we found B oc M' with 
&> 0.5, higher energy electrons should be accelerated in 
smaller clusters to emit synchrotron radiation at a given 
frequency. 

Deep radio survey with future radio telescopes (LOFAR, 
LWA, SKA) are required to test the presence of this cut- 
off/flattening in the luminosity function of the GRHs. 

The predicted evolution of the RHLFs with redshift is 
illustrated in Fig. 1131 the comoving number density of 
GRHs decreases with redshift due to the evolutions of the 
cluster mass function and of the probability to form GRHs. 
The decrease with redshift of the RHLFs calculated by 
adopting sublinear scaling of the magnetic field with cluster 
mass is faster than that in the superlinear scaling causing a 
spread in the RHLFs bundle with z. 

• Number counts (RHNCs) 

In Sec. 7 we have derived the integral number counts of 
GRHs at 1.4 GHz. We find that the number counts predicted 
for the same set of RHLFs discussed in Sees. 6 generally agree 
with those derived from the NVSS at the limit of this survey 
and within z — 0.2 fFig ll4|l . The fiattening of the counts 
below ~ 50 — 60 mjy is both due to the combination of the 
low power cut-offs of the RHLFs with the redshift limit, and 
to the RHLFs evolution with redshift. On the other hand, 
past calculations which assume a fixed fraction of GRHs 
with cluster mass predict an increasing number of sources 
at lower fiuxes (e.g., Enfilin & Rottgering, 2002). 

GRHs around the peak of our LFs {Piaghz ~ 10^*W/Hz) 
and at z~0.15 would be detectable at fiuxes below about 
20 mJy, which however is below the sensitivity limit of the 
NVSS for this type of objects. We estimate that the number 
of GRHs below this fiux could be up to 30-40 (whole sky, 
z ^ 0.2) if superlinear scalings of the mass with B hold. 

The predicted number of GRHs (Fie llSH (whole Universe) 
could be up to ^100 if a superlinear scaling of the mass with 
B holds, while a sublinear scaling would give a number 2-3 
times smaller. A substantial number of these objects would 
be found also down to a fiux of a few mJy at 1.4 GHz in the 
case of a superlinear scaling, while in the case of sublinear 
scalings the number of GRHs below about 10 mJy would be 
negligible. 

We also find that the bulk of GRHs is expected at 
z ~0.1-0.3 ('Fig |16ll . It should be mainly composed by those 
RHs populating the peak of the RHLFs, i.e. objects similar 
(or slightly more powerful) to the GRH in the Coma cluster. 

• Toward expectations at low radio frequencies 

In Sec. 7 we have extended our estimates to the case of 
low frequency observations which will be made with upcom- 



ing instruments, such as LOFAR and LWA. Lower energetic 
electrons contribute to these frequencies and thus - in the 
framework of the particle re-acceleration scenario - the effi- 
ciency of producing GRHs in galaxy clusters is expected to 
be higher than that of GRHs emitting at 1.4 GHz. 

By presenting the analysis for a representative set of pa- 
rameters, we have shown that the probability to have GRHs 
emitting at 150 MHz is significantly larger than that of 
those emitting at 1.4 GHz, particularly in the mass range 
~ 5 • 10^^ — 1.5 ■ 10^^ A/0. Consequently, the low mass cut- 
off in the RHMFs is shifted down by a factor of ~ 2. This 
is naturally expected and is due to the fact that slightly 
less turbulent systems are able to generate GRHs at lower 
frequencies. 

We have also estimated that the number counts of GRHs 
at low frequencies might outnumber those at 1.4 GHz by 
at least one order of magnitude. We venture to predict that 
LOFAR is likely to discover >10^ (all sky) GRHs down to a 
ffux of few mJy at 150 MHz. 
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